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Abstract — The work of Avestimehr et al. '07 has recently 
proposed a deterministic model for wireless networks and 
characterized the unicast capacity C of such networks as the 
minimum rank of the adjacency matrices describing all possible 
source-destination cuts. Amaudruz & Fragouli first proposed a 
polynomial-time algorithm for finding the unicast capacity of 
a linear deterministic wireless network in their 2009 paper. In 
this work, we improve upon Amaudruz & Fragouli's work and 
further reduce the computational complexity of the algorithm by 
fully exploring the useful combinatorial features intrinsic in the 
problem. Our improvement applies generally with any size of 
finite fields associated with the channel model. Comparing with 
other algorithms on solving the same problem, our improved 
algorithm is very competitive in terms of complexity. 



I. Introduction 

The deterministic channel model for wireless networks 
proposed by Avestimehr, Diggavi and Tse HI Q (referred 
to as ADT model thereafter) has been a useful tool for under- 
standing the fundamental limitations of information transfer 
in wireless networks. The ADT model captures two main 
features, the broadcasting and interference, that are present 
in wireless networks. It converts the wireless networks into 
deterministic networks, by making appropriate assumptions, 
that in turn lead to approximate capacity results. 

Consider a point-to-point Gaussian channel given by y = 
VSNRx + z where z ~ A/"(0,1) {M represents Gaussian 
distribution). Assume x and z are real numbers, then we can 
write y « 2" YJU *(*)2 _< + E£iO=(* + n) + z{i))2~* 
where n = [^logSNR] (here we assume a peak power of 
1 for x and z). If we think of the transmitted signal x as 
a sequence of bits at different signal levels, then the ADT 
model truncates x and passes only its bits above noise level 
(the first n most significant bits here), i.e., it converts the 
original Gaussian channel into a deterministic channel without 
noise. When applying the ADT model to wireless networks, 
the broadcasting is captured by the fact that in the resultant 
deterministic networks, all outgoing edges from the same 
signal level of any transmitting node carry the same unit 
information, and the interference is captured by the fact that 
at each signal level of any receiving node, only the modulo 
sum of all the received signals is available to the receiving 
node. This model is called the linear finite-field deterministic 
channel model in 0] iff). We refer to it as the ADT model 
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and denote the finite field of size p associated with the ADT 
model as F p in this paper. 

In fT) |2), the unicast (i.e., with one source S and one 
destination D) capacity C of any linear deterministic wireless 
relay network was characterized as the minimum rank of the 
adjacency matrices describing all its S-D cuts. An exhaustive 
search for finding the minimum rank of the adjacency matrix 
for all S-D cuts results in an algorithm with complexity 
exponential in the size of the network. 

Amaudruz & Fragouli [0 were the first to propose a 
polynomial-time algorithm for finding the unicast capacity of 
a linear deterministic wireless relay network (see also J4)). In 
this work, we improve upon Amaudruz & Fragouli's work and 
further reduce the computational complexity of the algorithm 
by fully exploring the useful combinatorial features intrinsic in 
the problem. Our improvement applies generally with any size 
of finite fields ¥ p associated with the ADT model. Comparing 
with other algorithms on solving the same problem f5) J6), our 
improved algorithm is very competitive in terms of complexity. 

This paper is organized as follows. In Section UT1 we briefly 
introduce the polynomial-time algorithm by Amaudruz & 
Fragouli for finding the unicast capacity of linear deterministic 
wireless relay networks. Section|III]gives a detailed description 
of our improvement upon the algorithm. First we introduce 
our improvement with an emphasis on the new components 
of our algorithm and how they fix the problems within the 
original algorithm. Then we explore several useful combina- 
torial features intrinsic in the problem. Finally we explain how 
these combinatorial features can be combined with our new 
components to reduce the complexity of the algorithm. We also 
give the comparison results between our improved algorithm 
and other algorithms on solving the same problem. Section ITVl 
concludes the paper. 

II. Preliminaries and Background 
A. Notations and Definitions 

In 12, it is shown that an arbitrary deterministic relay net- 
work can be expanded over time to generate an asymptotically 
equivalent (in terms of transmission rate) layered network. 
Therefore, we focus on layered deterministic networks. 

Let Q = (V, £) denote a layered deterministic wireless relay 
network where V represents the set of nodes in the original 
wireless relay network, each node in V has several different 
levels of inputs and outputs and £ is the set of directed edges 
going from one input of some node to one output of some 



other node. For example, Fig. |l(a)| gives a graph representation 
of a layered deterministic wireless relay network where each 
node is labeled with a capital letter, all inputs (outputs) from 
nodes are labeled as {x^ ({yj}), 1 ^ i, j ^ 8. In the layered 
network Q, all paths from the source node S to the destination 
node D have equal lengths J2] . The set of nodes V are divided 
into different layers according to their distances to S. The first 
layer consists of S and the last layer consists of D. Let A{xi) 
(or A(yj)) denote the node where an input Xi (or an output 
Dj) belongs to. Let C(A) (or C{xi), £{yj)) denote the layer 
number where node A (or Xi, yj) belongs to. Denote M as the 
maximum number of nodes in each layer, L the total number 
of layers and d the maximum number of outgoing edges from 
any input in any node in the network Q in this paper. 

A cut ft in Q is a partition of the nodes V into two disjoint 
sets SI and il c such that S £ il and D £ il c . A cut is called a 
layer cut if all edges across the cut are emanating from nodes 
from the same layer, otherwise it is called a cross-layer cut. 
An edge (xj, yj) £ £ belongs to layer cut I if C(xi) = I. 

The adjacency matrix T(x,y) for the sets of inputs x = 
{x 1 ,x 2 ,...x m } and of outputs y = {y 1 ,y 2 , ...y n } in Q is a 
matrix of size m x n with binary {0, 1} entries. The rows 
correspond to {xi £ x} and columns corresponding to {yi £ 
y} and T(i,j) = 1 if (xi,yj) £ £. The adjacency matrix 
T(E) for a set of edges, E, is the adjacency matrix for the 
sets of their inputs and their outputs. 

A set of edges, E, are said to be linearly independent (LI) if 
rank(T(_E)) = (where the rank is computed over GF(2)), 
otherwise they are said to be linearly dependent (LD). In Q, 
each S-D path is of length L — 1 and crosses each layer cut 
exactly once. A set of S-D paths are said to be LI if the subsets 
of their edges crossing each layer cut are LI, otherwise they 
are said to be LD. In this work, we will consider a slightly 
more general adjacency matrix, where the non-zero entries can 
be from a finite field T p , and the rank is also computed over 
J- p . Of course, all our results will also apply to the binary 
field case. 

Let £q be the set of edges crossing the cut f2 in Q. The 
cut value of is defined as rank(T(£fj)), which based on 
the definition equals the maximum number of LI edges in £q. 
Note that the cut value defined above is different than that for 
regular graphs (which is just the number of edges crossing the 
cut). It is proved flTJ fl2) that the unicast capacity of a linear 
deterministic wireless relay network is equal to the minimum 
cut value among all S-D cuts. 

B. Algorithm by Amaudruz & Fragouli 

The unicast algorithm by Amaudruz and Fragouli finds 
the maximum number C of linearly independent S-D paths in 
a given layered linear deterministic relay network Q, where 
C is the unicast capacity of the network. The algorithm is a 
path augmentation algorithm, operating in iterations. In each 
iteration, the algorithm tries to find an additional S-D path 
so that all S-D paths found are LI. Let V = {Vi,...,V k } 
denote the set of k LI S-D paths found in the first k iterations. 
In the process of finding the (k + l)-th S-D path Vh+i in 
iteration k + 1, the algorithm may make modifications to V 



while still maintaining a set of k LI complete S-D paths. The 
unicast algorithm determines Vk+i by exploring nodes in Q 
in a certain order as outlined shortly. 

The algorithm is implemented in two recursive functions 
Ea and E x that explore a node and input respectively. The 
exploration of a node A takes place when Vk+i has been 
extended from S to A and needs to be completed from A to 
D. In iteration k + 1, the unicast algorithm calls Ea with the 
following inputs: Q, V = {Pi, ...,Vk}, the indicator function 
Ai (that implements a marking mechanism for visiting nodes 
and inputs/outputs) and S. The function Ea returns true with 
one more S-D path Vk+i recorded in V if it succeeds in finding 
Vk+i, false otherwise. 

Exploring node A implies exploring all unused inputs {xi} 
of A. So we explain the exploration of an input Xi of A below. 
Hereafter, denote U as the sets of used edges by V in layer 
cut I and XJ X and U l y as the sets of inputs and outputs used 
by U l . Let C(xi) = I. If Xi £ U l x , do nothing. Otherwise, 
consider each yj with (xi,yj) £ £ as follows. 

(a) yj is used. Let L Xi denote the smallest subset of U l x with 
s = \L Xi \ ^ \U l x \ = k such that T({L Xi ,Xi}, U l y ) has 
rank s. The authors prove that replacing any Xk £ L Xi 
with Xi, the algorithm can still maintain k LI S-D paths 
and the task now is to complete Vk+i from A(xk)- So in 
this case the unicast algorithm first finds the set L Xi in 
function FindL. Then it replaces each Xk £ L Xi with xi 
and calls a Match function to find a new set of k edges 
in layer cut I to maintain k LI S-D paths in V and tries 
to complete Vk+i from A(xk) if A(xk) is not marked 
or from Xk if Xk is not marked. We refer to this step as 
same-layer rewiring. 

(b) yj is not used. A rank computation function is called on 
the matrix T({U X , Xi}, {U y , yj}). If the matrix is not full 
rank or A(yj) has been visited before, do nothing. If the 
matrix is full rank and A(yj) has not been visited before, 
add (xi,yj) to Vk+i and try to complete it from A(yj) by 
exploring A(yj). We refer to this step as forward move. 
If it fails to complete Vk+i from A(yj), a 0-function is 
called for each yk £ U y with A(yk) = A(yj). Let V Vk be 
the path using yk and let (xk,yk) £ U l be the path edge. 
The idea of the 0-function is to complete Vk+i from 
A(yj) to D using the partial path of V yh from A(yj) to D 
and then try to complete the path V yk from A(xk). The 0- 
function does the following: remove {xk,yk) from the set 
of used edges and try to complete V Vk from A(xk). We 
refer to this step as backward rewiring. The 0-function 
will be executed at most M times. 

We refer the reader to O for more details. The complexity 
of the algorithm is 0(M -\£\ - C 5 ) and its computational parts 
include the FindL, Match and rank computation functions each 
with complexity 0(fc 4 ), 0(k 3 ) and 0(k 3 ) respectively. 

C. Other Related Algorithms 

Yazdi & Savari |5| developed another polynomial time 
algorithm with complexity 0(L 8 M 12 hl + LM 6 CIiq) (where 
ho denotes the maximum total number of inputs/outputs at any 
layer) by relating matroids with this problem. Most recently, 
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Goemans, Iwata and Zenklusen (6) proposed a strongly poly- 
nomial time algorithm for this problem, whose complexity is 
0(LM 3 log M), i.e., it does not depend upon C. 

III. Improved Unicast Algorithm 

In this section we outline certain improvements that can be 
made to the algorithm of J3). In particular, we elaborate on 
several useful combinatorial aspects that allow us to reduce the 
overall time complexity. Moreover, these improvements also 
fix certain issues with the original algorithm 0. As mentioned 
previously, our proposed improvements apply over arbitrary 
finite fields. 

A. Improving the Original Algorithm 

The main idea in Q is to find path Vk+i in iteration k + 1 
while maintaining linear independence among all S-D paths in 
V. In this process, previous paths may be rewired. However, 
there are cases when the original algorithm may fail to find the 
exact unicast capacity. We illustrate this using the following 
examples. We point out that these issues seem to have been 
resolved in HJ. However, our proposed algorithm has several 
differences from J?) as discussed at the end of Section UlI-DI 
Improved Backward Rewiring 

We use the example in Fig. [TJ to show that there are cases 
where the 0-function above is insufficient, causing failures of 
the original algorithm. Then we illustrate how it can be fixed 
by introducing an improved backward rewiring mechanism. 

In Fig. | l(a)| three LI S-D paths with color red, green and 
blue are found in the first three iterations of the algorithm. 
Let's see how the algorithm goes in iteration four. Let's say 
the algorithm has extended Va along the purple path to 2/20- 
The call E A (G, V,M,N) fails since the only input a; 24 of N 
is used by paths in V . So 0-function is called on y%g and then 
node I is explored in Ea(G, V,Ad,I), but since there is only 
one path from all inputs of I to D, Ea{G,V,AA, I) fails, and 
finally the algorithm returns false and reports unicast capacity 
of 3. However, the unicast capacity of the network is 4 and a 
capacity-achieving transmission scheme is given by the four 
S-D paths in Fig. 1 1 (b)| in different colors. 

We propose the following improved backward rewiring 
mechanism to fix the problem above and to replace the original 
0-function. Let A denote a node in the network (not to be 
confused with A in the figure). First, the backward rewiring 
is allowed on every node A whenever it is explored in finding 
Vk+i- Second, the backward rewiring on node A includes the 
following operations. Let C(A) =1+1, For any output y 
of A with y £ Uy and y is used by a path in V at the 
beginning of the current iteration (if such y exists), (1) find 
one x € U l x such that T(U X - x, XJ l y - y) has full rank, (2) 
then rematch (U x — X,Uy — y) to generate a new set of k LI 
used path edges in layer cut I and (3) finally try to complete 
the partial path from A(x). Lemma [3] guarantees that for a 
given y 6 U l y there is always one such x and also a set of 
edgefO] P y _> x = = y), (2:1,2/2), (£2,2/2), (2:2,2/3), 

'We use the notation P y ^ x since this set of edges can be interpreted as 
an alternating path, as we show in Section Ull-B I 



...(aw_i,j/ m /), {x m i = x,y m >)} = {e x , e 2 , e 2m .'_i} with 
< i < m' being edges used by V, which can 
be found with complexity 0(k 3 ) and 0(k 2 ) respectively. 
Along the alternating path P y ^ x , the rematching of the used 
path edges in layer cut I can be done easily as follows: 
U l = U l - ei + e 2 - e 3 + ... - e 2m '-i- 

Consider applying our improved backward rewiring in the 
example in Fig. [TJ It happens on the outputs of nodes N and I. 
Its application to N is straightforward. Let's look at its appli- 
cation at the output 2/14 of node I. First it finds xq <E U 2 with 
T(U 2 —xe,Uy— 2/14) having full rank and the alternating path 
p yi4->x 6 = {(x 7 , 2/14), (257 1 2/13), (xe, 2/13)}- The rematching is 
done by U 2 = U 2 - (x 7 ,y u ) + (2:7,2/13) - (^6,2/13)- Then 
node B = A{xq) is explored. Finally the improved algorithm 
returns four LI S-D paths in Fig. |l(b)| as expected. 
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Fig. 1. Illustrating example for improved backward rewiring 

Improved Same-Layer Rewiring 

We use the example in Fig. [2] to show that the same-layer 
rewiring in original algorithm is insufficient. Suppose the red 
S-D path is found in the first iteration. In iteration two, suppose 
that the algorithm first extends V2 along the green path to x 4 . 
The same-layer rewiring from 2:4 will mark X3. Since T(x$ + 
2:4,2/5 +2/6) is not full rank, the algorithm fails to complete 
■p 2 along the green path. It continues to extend V2 along the 
blue path to x&. Since 2:3 is marked, the same-layer rewiring 
from 2:5 won't be applied on X3 and the call Ea(G, P,A4,C) 
fails. The algorithm finally returns false and reports unicast 
capacity of 1, However, the network has a unicast capacity of 
2 indicated by the two paths in Fig. |2(b)| 

We develop our improved same-layer rewiring to fix the 
above problem as follows. First, an input Xk should not 
be blocked from being visited via same-layer rewiring from 
any input x. t just because it has been visited via same-layer 
rewiring from another input Xj, Consider the example in Fig. 
12 If we allow 2:3 to be visited via same-layer rewiring from X5, 
the algorithm may succeed in finding two LI paths as indicated 
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in Fig. |2(b)| However, this needs to be done carefully. Consider 
again the example in Fig. 12 If we allow same-layer rewirings 
from all inputs, then we might run into an infinite loop of 
going from x$ to x% via same-layer rewiring and going from 
X3 to £5 via same-layer rewiring and so on. 

The goal of a same-layer rewiring operation in iteration k+1 
is to ensure that every input, which allows the algorithm to 
maintain k LI S-D paths and can further extend the current 
partial path, has the opportunity of being explored, while 
ensuring that we do not enter an infinite loop. In this work 
we achieve this by using a pair of labels of each node. 




(a) (b) 



Fig. 2. Illustrating example for improved same-layer rewiring 

Each node has a label that takes values - "explored" or 
"unexplored". The other label is a type that takes values 1, 2. 
We initialize the type of every node to be 1 at the beginning of 
the iteration. A type 1 input is allowed to initiate same-layer 
rewirings. An input that is explored via a same-layer rewiring 
from a type 1 input x t is assigned as type 2. A type 2 input 
is not allowed to initiate same-layer rewirings to avoid the 
possibility of infinite loop. If an input x (of either type) is 
explored via a backward rewiring, it is re-assigned as type 1 
(since U l x and U y change since last time x was explored). 

Consider applying our improved same-layer rewiring in the 
example in Fig. [2] X3 is first visited via a same-layer rewiring 
from X4 (of type 1) when it is assigned as type 2. Later on 
X3 is revisited via a same-layer rewiring from x§ (of type 
1) when it is assigned as type 2 again, so it won't initiate a 
same-layer rewiring to 2:5, instead it only looks for a possible 
forward move which happens along the edge (#3, 2/5) (and the 
improved algorithm finally succeeds in finding 2 LI paths as 
in Fig. 

B. Useful Combinatorial Features 

In this subsection, several useful combinatorial features 
intrinsic in the problem are introduced which are used later in 
our improved algorithm to reduce the complexity. 

In the following, we define a set A Xj similar to but more 
general than L Xi in the original algorithm by Amaudruz and 
Fragouli. A Xi applies to any size of finite field F p associated 
with the ADT model for the network. 

Definition 1: Define A Xi as a subset of U X when Xi is 
explored such that 

T(x it UfM) = J2 <-T( Xj ,U y c ^). (1) 

where {a Xi } are non-zero coefficients from F p . 

Lemma 1: A Xi and the set {a Xi } are unique and can be 
found with complexity 0(k 3 ) in iteration k + 1. 



Since T(U X (X '> ,Uy (x ' } ) has full-rank, A Xi and the set 
{a Xi } are unique and can be found with complexity 0(k 3 ) 
by using Gaussian elimination. 

C(x ) 

Let y x . denote the bipartite graph containing nodes U x U 
U y when Xi is explored in iteration k + 1 and Q x . denote 
the bipartite graph containing nodes {x^ U U x Xt ' U Uy*'. 

In the following, we refer to an alternating path as a path in 
which the edges belong alternatively to the set of used edges 
and the set of unused edges. 

Lemma 2: There is an alternating path from Xi to any 
Xj G A Xi in the graph Q+. of the form P Xi ^ Xj = {(xuyi), 
(xi,yi), (xi,y 2 ), (2:2,2/2), ..-(im-hJ/m), (x m = Xj,y m )} 
with (x g ,y q ),l < q < m being edges used by V. The 
complexity for finding these \A Xi | paths is bounded by 0(k 2 ) 
in iteration k + 1. 

Proof: Let C{x l ) = I. Given rank(T(i7*, U l y )) = k, 
for any x 3 G A Xi , rank(T '(U l x ,U y )) = rwak(T{U l x + x % - 
Xj,U l y )) = k where k = \V\ in iteration k + 1. Introduce 
an auxiliary output y' and an edge (xj,y'). It's easy to see 
that rank(T([4 + x { , U y + y')) = k + 1. Let G++ denote the 
bipartite graph containing nodes {.t,} U U l x U U y U {y 1 }. 

Given T(U X , U y ) has full rank, we know that the polynomial 
of the determinant of the Edmonds matrix of the bipartite 
graph Q x . is not identically zero, so there is a size k perfect 
matching in Q Xi 0, Mi = U l giving such a matching. 
Similarly given Twk(T(U x + Xi,U y + y')) = k + 1, there is a 
size k+1 perfect matching in G Xi + - By Berge's Lemma (8), we 
know that there is an alternating path, relative to the matching 
Mi, starting from an unused input Xi to an unused output 
y', alternating between edges not in the current matching Mi 
and edges in the current matching Mi, i.e., there is a path 
Px^y' = {(x l ,y 1 ), (xi,yi), (x 1 ,y 2 ), (2:2,2/2), ...(x m -i,y m ), 
(x m ,y m ), {x m = Xj,]/)} with (x q ,y q ),l < q < m 
being edges in Mi. So we proved that there is an alter- 
nating path P Xi ^ Xj = {(xi,yi), (xi,y%), (2:1,2/2), (2:2,2/2), 
■~{x m -i,y m ), (x rn = Xj,y m )} with (x q ,y q ),l < q < m 
being edges in Mi = U . 

Since the number of nodes in (/+ is bounded by O(k), the 
number of its edges is bounded by 0{k 2 ). Finding P Xi ^ x for 
all Xj G A Xi in can be done with complexity 0(k 2 ) with 
some well-known graph traversal algorithms, like breadth-first 
search J9). ■ 

Lemma 3: Let rank(T(U l x ,U y )) = \U l x \ = \U l y \ = k + 1. 
Given any y G U y , there exists at least one x G U x , such that 
mnk(T(U x — x,U y — y)) = k. Moreover there is an alternating 
path from y to x of the form P y ^ x = {(xi,yi = y), 
(2:1,2/2), (2:2,2/2), (2:2,2/3), -(x m '-i,y m >), (aw = x,y m >)} 
with (x q , y q ),l < q < m! being edges in U l . The complexity 
of finding one such x is bounded by 0(k 3 ) and the complexity 
of finding path P y ^ x is bounded by 0(k 2 ). 

Due to lack of space, we skip the proof here. The proof 
of existence of P y ^ x is similar to Lemma |2] by introducing 
an auxiliary input x' and output y' and edges (x',y), (x,y') 
leading to rank(T(L^ + x\ XJ\ + y')) = k + 2. 

Lemma|4]develops an equivalent but computationally simple 
method to speed up the rank computation when x^ is explored 
given A Xi and the set of associated coefficients {a Xi }. 
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Lemma 4: Let T(U x ,U y ) have full rank k. The rank com- 
putation for checking mnk(T(U x + Xi, U y + y)) = k or k + 1 
for any x 4 ^ {/J., £(2;*) = Z, y ^ {/' and £ S is 

equivalent to checking T(xi,y) = Ex-gA • T(xj,y) or 
not, with complexity bounded by 0(k) given A Xj and {a Xi }. 

Proof: Given T(U l x ,U l y ) has full rank fc, rank(T(?7j. + 
Xi,U y + y)) = k is equivalent to that T{xi 1 U l y + y) = 
Sx^eAj,. fl i; ' T ( x j, u i + V) for some A i, £ and {a x ,}. 
Since A Xj C U l x and the set {ajj are unique for which 

T{x l ,U l y ) = ExgA^ < ■ T ( x j> ul y) holds ( b y Lemma [D, 
there must be A' x . = A Xi and {a Xi } — {a x >}. This leads 
to that rank(T(U l x + Xi 7 U y + y)) = k is equivalent to 
T(xi,y) = E^gA.. < ' T ( x hV)- ■ 
Lemma 5: Let x' € A Xi . If x' is explored via a same- 
layer rewiring from x%, A x i = A Xi + Xi — x' and the set 
of associated coefficients {a x /} can be computed from {a Xi } 
with complexity 0(k) in iteration k + 1. 

Proof: Let £(x.j) = I. Note that when x' is explored via 
a same-layer rewiring from Xj, U l x is updated as U l x — x' + Xi, 
U l y is unchanged and T(U l x — x' + Xi, U l y ) has full rank. Based 
on definition, 

T{ Xl ,U y )= Y, ^■T{x^U l y )+a' x ^T{x',U l y ). (2) 

a:jGA x .\a;' 

where {a Xi } are non-zero coefficients from W p . So we have 

T ( x '^ y ) = J2 ^ ± - T ( x j,U l y )-^-T(x i ,U l y ). (3) 

Xj£A Xi \x' x * Xi 

Since T(U l x — x' + Xi, U l y ) has full rank, equation (O is the 
unique way that the row T(x', U l y ) can be expressed as a linear 
combination of the rows in this matrix. So we conclude A x i = 
A Xi + Xi — x' and the set of associated coefficients {a x > } can 
be computed from {a Xi } with complexity O(k). Note that in 
iteration k + 1, \A Xi | < \U l x \ = k. ■ 

C. Reducing the Complexity and the Overall Algorithm 

As mentioned before, the computational parts of algorithm 
J5) include the FindL (finding L Xi ), Match (update U after a 
same-layer rewiring from Xi) and rank computation functions. 
Now we explain how the combinatorial features from Section 
IIII-BI can be used to further reduce the complexity of the 
unicast algorithm. 

Lemma Q] shows that A Xi and the set of associated coeffi- 
cients {a Xi } for any type 1 input Xj can be computed with 
complexity 0(k 3 ) in iteration k + 1. Lemma [5] tells that for 
any type 2 input x', x' € A Xi , that is explored via a same- 
layer rewiring from a type 1 input xu A x / and the set of 
associated coefficients {a x '} can be computed with complexity 
O(k) given A Xi and the set of associated coefficients {a Xi }. 

Second, based on Lemma [2] the matching or updating of 
U after same-layer rewirings from any type 1 input Xi can 
be done with complexity 0(k 2 ) in iteration k + 1 as follows. 
First find all [ A^,. j paths P Xi ^ Xj , Vxj £ A Xi with complexity 
0(k 2 ) for Let P Xi ^ X] = {(xi,yi), (xi,yi), (x 1 ,y 2 ), 
—(x m -i,y m ), ( x m = x j,y m )} = {ei,e 2 ,...,e 2m } with 
(x q , y q ),l < q < m being edges used by V for any Xj G A Xi . 



Then updating of U^ Xi > after a same-layer rewiring from Xj 
to Xj can be done by U c ^ Xi ^ <— U ctyXi ^ + e\ — e 2 + ... — e 2m . 

Third, Lemma|4]tells that the rank computation in a forward 
move from any x- t (either of type 1 or of type 2), Xi ^ U l x , 
C(xi) =1, for checking Twk(T(U x +Xi,U y +y)) = k or k+1 
for any y U y and (xi,y) e £ is equivalent to checking 
T(xi,y) = Ex eA,. a i, ' T ( x j^v) or not ' with complexity 
bounded by O(k) given A Xi and {a Xi } in iteration k + 1. 

Finally, as mentioned before, in our improved backward 
rewiring from an output y, to find one x with T(U l x —x, U l y —y) 
having full rank and to rematch (U x — x,U y — y) can be done 
with complexity 0(k 3 ) in iteration k+1 guaranteed by Lemma 
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Table U gives an overall description of our improved 
unicast algorithm which is implemented in a function 
Ea(G,V,M.,A) where all inputs are the same as in the 
original algorithm. A complete software implementation of our 
improved unicast algorithm can be found in iflOl . 

TABLE I 

Pseudo-code for our improved algorithm 



{(T,F)\=E A (G,V,M,A) 
' M(A) = T,C(A) = I 

U l = {used edges in layer cut I}, U' x = {x t € U 1 }, U' y = {y :j G U 1 } 
for any x : A(x) = A. x U l x , M(x) = F, GetType(a;) = 2 
' M(x) = T 

for any y : (x, y) G £,y $ Uy, M(A{y)) — F //forward move 

■ ifT(x,y)^-£ X3eAx ai-T(x j ,y) 

{ f Update (P) ; £/' <- U l + e 
< I if A(y) = D, return (T) 

1 else if E A (Q, V, M, A(y)) = T, retum(T) 
{ [ U' <- U l - e; Restore(P) 
for any x : A(x) = A, x U' x , M(x) = F, GetType(a;) = 1 
' M{x) = T 

Compute A x and the set of coefficients {a x } 

for any y : (x, y) G £,y Ul, M(A(y)) — F //forward move 

■ XT(x,y) ^-£^^4 -T( X: j,y) 

IUpdateCP); U l <- U l + e 
if A(y) = D, return (T) 
else if E A (Q, V, M, A(y)) = T, return(T) 
U l <~U l - e; Restore^) 
' < Find all paths P x ^ Xj for all \/ Xj G A x 

for any Xj : Xj e A x with P x ^ Xj = {ei, e 2 , ...e 2m } = 

{(x, yi),(xi,yi), (2:1,2/2), ---(xm = Xj,y m )} //same-layer rewiring 

■ M(xj) = F; SetType(xj, 2); 
A Xj = A x - Xj + x 

compute {a Xj } based on {a x } according to Lemma [5] 
] Update(P); U l <- U l + ei - e 2 + ... + e 2 ,„-l - e 2m 

if E A (g,V,M,A(xj)) = T, return(T) 
I U l <- U l - ei + e 2 - ... - e 2m _i + e 2m ; Restore(P) 
for any y : A(y) = A,y G U^Miv) = F 

and y is used by V at the beginning of the iteration //backward rewiring 
( M(y) = T 

find one x G C4 _1 with TiU^r 1 - x, U^ 1 - y) having full rank 

and find P v ^ x = {ei,e 2 , ...e 2m /_ 1 } 
J = {(xi,yi = y), (xi,y2), (2:2,3/2), ■■•(a: m ' = a:,y m /)} 
] M{x)= F, SetType(K, 1) 

Update(P); U 1 ' 1 <- C/'- 1 - ei + ej - ... - e 2m <_i 

If V, M, A(x)) = T, return (T) 

. U'- 1 <- + ei - e 2 + ... + e 2m /_ 1 ; Restore(P) 

return (F) 



D. Complexity Analysis and Comparison with Existing Results 

To analyze the complexity, we first bound the total number 
of inputs of different types being visited in each iteration 
of the algorithm. Note that once a node or input/output is 
visited/explored, it's labeled as explored (by Ai) and not 
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allowed to be explored again unless it is relabeled as unex- 
plored again. At the beginning of each iteration, all inputs 
are initialized as unexplored type 1 inputs whose number is 
bounded by 0(|V X |) (let V x ={all inputs in the network}). In 
each backward rewiring operation, one input will be assigned 
as unexplored type 1 input. From the definition of backward 
rewiring, the total number of valid outputs that initiate a 
backward rewiring is no more than | Vcc |, which means the 
total number of backward rewiring operations is bounded by 
0(| V x I)- So the total number of type 1 inputs being visited 
is bounded by 0(|V X |) in each iteration. In each same-layer 
rewiring operation from a type 1 input, one input will be 
assigned as unexplored type 2 input. The total number of same- 
layer rewiring operations from any type 1 input x is no more 
than \A X \ ^ k in iteration k+1. So the total number of type 2 
inputs being visited is bounded by 0(A;|Va;|) in iteration k + 1. 

The worst case in computation in iteration k+1 are no more 
than: (1) for each type 1 input Xi, compute A Xi and {a Xi } with 
complexity 0(fc 3 ) and find all paths P Xi ^ x for Vxj G A Xi 
with complexity 0(k 2 ), (2) for each type 2 input Xj, compute 
A Xj and {a Xj } with complexity O(k), (3) for each type 1 
or type 2 input x, compute rank for T{U l x + x, U l y + y) for 
all y U l y , (x, y) E £ with complexity 0(k) given A x and 
{a x } (for any x, the total number of such y is no larger than 
d) and (4) in each backward rewiring from a certain y, find 
one x with T(U l x — x,U y — y) having full rank and to rematch 
{U x — x, U l y -y) with complexity 0(fc 3 ). Note that k < C. It's 
obvious that the total complexity of our improved algorithm 
is bounded by 0(114 1 • C* 4 + d ■ \V X \ ■ C* 3 ). 

Due to lack of space, we skip the proof of correctness 
for our improved algorithm, however a complete and detailed 
proof can be found in iflOl . 

Table HI1 lists the comparison results between different algo- 
rithms for finding the unicast capacity of linear deterministic 
wireless relay networks, specially in their complexity. 

TABLE II 

Comparison of algorithm complexity 



Algorithm 


Complexity* 


Notes 


|3| 

El 


0(M\£ |C D ) 
0(d\V*\C 5 + \V y \C 5 ) 


Always higher than ours 
especially when C is large 




0{L a M rz h*+LM Q Ch%) 


Always higher than ours, especially when M or L 
is large 


m 


OiL^M*- log(ML)) or 
0(LM 3 log M) 


Straightforward comparison is not possible. 6j will 
have lower complexity if C is much larger than M 


Our work 







* Denote C as the unicast capacity, M the maximum number of nodes in each layer, L the total 
number of layers, d the maximum number of inputs of any node, ho the maximum number of 
inputs/outputs at any layer, E the total number of edges, V* the total number of inputs and 
V y | the total number of outputs. Note that M > d (since by definition each input can have at 
most one connection to each node in the next layer), £ > (because of broadcasting) and 
ho > C (based on definition). 

We note that the issues with the original algorithm Q 
mentioned in Section IHI-AI have been fixed in @]. The main 
difference between our improved algorithm and the algorithm 
in J4) is that our improved algorithm utilizes those useful 
combinatorial features intrinsic in the problem described in 
Section IIII-BI which lead to reduced complexity. The other 
difference comes from the same-layer rewiring and backward 
rewiring. In (4), the same-layer rewiring starts on each input 
at most once (using the ML indicator function) while our 
algorithm allows multiple same-layer rewirings starting from 



certain inputs (that is, if an input is explored via a backward 
rewiring, it is reassigned as type 1 input and allows to initiate 
same-layer rewiring again). In [4], the backward rewiring 
(implemented in (^-function there) allows exploration on every 
Xfc G U x such that the resulting adjacency matrix of used path 
edges still remains full rank while our algorithm only finds 
one such Xk G U x and explores it. Note that it can be verified 
that the combined effects of the different same-layer rewiring 
and backward rewiring in two algorithms are the same. 

IV. Conclusions 

An improved algorithm for finding the unicast capacity 
of linear deterministic wireless networks is presented. Our 
algorithm improves upon the original algorithm by Amaudruz 
& Fragouli. We amend the original algorithm so that it finds 
the unicast capacity correctly for any given deterministic net- 
works. Moreover we fully explore several useful combinatorial 
features intrinsic in the problem which lead to reduced com- 
plexity. Our improved algorithm applies with any size of finite 
field associated with the ADT model defining the network. 
Our improved algorithm proves to be very competitive when 
comparing with other algorithms on solving the same problem 
in terms of complexity. 
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